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We offer a revised exposition of the three types of heat-propagation theo- 
ries proposed by Green and Naghdi, in the form they were given in 1992-3 
[9] [10]. Those theories, which make use of the notion of thermal displace- 
ment and allow for heat waves, are at variance with the standard Fourier 
theory; they have attracted considerable interest, and have been applied in 
a number of disparate physical circumstances, where heat propagation is 
coupled with elasticity, viscous flows, etc. [27]. However, their derivation is 
not exempt from criticisms, that we here detail, in hopes of opening the way 
to reconsideration of old applications and proposition of new ones. 
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1 Introduction 

The nonconventional heat conduction theory by Green and Naghdi [8] [9] [10] 
provides a general framework within which a wider range of thermal prob- 
lems than within the standard theory have been modeled. The theory is 
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subdivided into three different types of different generality: Type I, under 
special circumstances, encompasses the standard theory based on Fourier's 
constitutive prescription for the heat flux vector; Type II allows for propa- 
gation of thermal waves without internal energy dissipation; and Type III is 
meant to provide a framework to describe an even wider range of problems, 
many of which are given a reasoned exposition in a recent book [27]. In 
fact, Green and Naghdi developed their theories, whatever the type, within 
a thermodynamic framework that they regarded as robust enough to deal, in 
addition to heat propagation in a rigid conductor, with a number of coupled 
phenomenologies such as those covered by thermoelasticity [10] or by the 
theory of thermoviscous fluids [11^ [T2] . 

This paper is devoted to point out and amend the faulty arguments and 
moot points we found in Green-Naghdi theories, none of them excluded. 
Since the reservations we have are of a general thermodynamic nature, we 
confine our discussion to heat propagation in rigid conductors. In Section 2 
we present our understanding of the thermodynamic structure common to 
all three types of Green-Naghdi theories. This section serves as an indis- 
pensable background for our type-by-type scrutiny, to come in Sections 3 to 
5. Our main criticisms are listed in the conclusive Section 6. 



2 The common structure of Green-Naghdi theo- 



We begin by recapitulating the field equations used by Green and Naghdi, 
in the form those equations are given in [10], a few notational changes apart. 

2.1 Entropy balance 

At variance with the classical procedure of postulating balance of energy (the 
so-called First Law) and imbalance of entropy (the Second Law), Green and 
Naghdi base their approach on a statement of entropy balance, a procedure 
they first introduced in [6] (the consistency of the latter approach with the 
former is demonstrated in [7J). Their entropy balance reads: 



with 77, s and £, respectively, the entropy, external entropy supply and in- 
ternal entropy production, and with h the entropy influx vector (here and 
henceforth a superposed dot signifies differentiation with respect time t). 
They assume that the energy inflow (q, r), where q is the heat influx vector 



ries 



i] = —divh + s + ^ 
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and r is the external energy supply, is proportional to the entropy inflow 
(h, s) via the absolute temperature 9: 

q = 9h and 9s = r, with 9 > 0. (2) 

Remark. The generality of this classic assumption has been questioned by 
Miiller [221 ED [23] and others, e.g. Liu HQ QI1 Q21 E21 123 • To the best of 
our knowledge, Green and Naghdi simply ignored those reservations, which 
however were considered briefly in [5]. The issue was taken up in [2] by 
Bargmann and Steinmann, whose developments and results were straight- 
ened and completed by us in our paper pQ. Both in [2] and in [I], the 
Green-Naghdi theory of Type III is considered within a Miiller-Liu frame- 
work; among other things, it is shown that ([2]^ holds if both h and q depend 
isotropically on the state variables, whereas it does not hold if that depen- 
dence is only transversely isotropic. That isotropy guarantees ©i, whereas 
transverse isotropy does not, was first suggested by Liu [20], for elastic bod- 
ies studied in a Miiller-Liu framework. Recently, one of us has offered a 
simple argument to prove, in a standard thermodynamic setting, that in- 
flux proportionality holds (does not hold) if heat conduction is isotropic 
(transversely isotropic) [26]. 

2.2 Energy balance, free energy, and the reduced entropy 
balance 

Multiplication by 9 of Eq. ([T]) yields: 

0rj = -div (q) + 9~ 1 q-V9 + 9{s + £) = ~divq + 9" 1 q ■ V9 + r + 9$. 
With this, the energy balance 

i: = — divq + r, (3) 
and the notion of Helmhholtz free energy per unit mass 

i/j = e-9rj, (4) 
Green and Naghdi arrive at the following reduced entropy equation: 

ip + r ] 9 + 9- 1 q-V9 + 9t = 0; (5) 
here 9^ is the internal dissipation. 
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2.3 The thermal displacement 

Relation ([S]) is the common point of departure of the constitutive devel- 
opments in Green-Naghdi Type I, Type II, and Type III, theories, to be 
recapitulated here below. In all those theories, (0) plays the same role as 
the reduced entropy inequality 

tp + r]9 + 6- l q- V6»<0 (6) 

in the standard Coleman-Noll procedure [3]. The difference is, of course, 
that Green and Naghdi do not assume once and for all that the internal 
dissipation should never be negative. Another major difference with the 
standard treatment of heat conduction is that, in all their three types of 
theory, Green and Naghdi include among the state variables, either implic- 
itly or explicitly, the thermal displacement a, that is, by definition, a time 
primitive of an empirical temperature T: 

a := T; 

furthermore, they take the absolute temperature 9 to be an affine function 
of T. Hereafter, for simplicity, we assume that empirical and absolute tem- 
perature coincide: 

a := 6. (7) 

Remark. A short historical account of the fortunes of thermal displacement, 
a notion that was first introduced by Helmholtz [13] in 1884, is found in the 
appendix of |24j . While the formal role of this notion in postulating the basic 
balance laws of thermomechanics is clear |24j . a physical interpretation is still 
wanted, especially if consistent with the statistical-mechanics interpretation 
for temperature. 

2.4 State spaces 

For a state, we mean a list of independent state variables, whose instanta- 
neous values at a space point determine the current local state of a material 
body during a thermodynamic process of interest. Green and Naghdi based 
their approach on the introduction of three state spaces, one for each of 
their three types of heat conduction theory. Unfortunately, they were not 
consistent with their choice of state spaces: in their first paper [8], state 
spaces of Types II and III include the thermal displacement a; not so in 
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their follow-up papers [9] and [10] . We have not been able to find any com- 
ment of theirs about this important change in their views. In the following, 
we adopt the lists of state variables that are found in [8], namely, 

Si = {9, V9} = {a, Vd}, <Sn = {a, d, Va}, Sm = {a, a, Va, Vd} 

(needless to say, each subscript points to the type of theory the state it 
individuates is chosen for). Note that, in view of J7|), is included in all 
three notions of state. Note also that Type III is the most inclusive theory, 
while Types I and II cannot be ordered with respect to generality. 

2.5 Thermodynamic restrictions on response mappings 

Assume, as is customary, that all of the mappings for the fields in ([5]) that 
are in the need of constitutive specifications (namely, tp, r], q, and £) depend 
in principle on one and the same chosen list of state variables S: ip = ip(S), 



Let us now choose S in (|8|) to be Si. Then, granted the arbitrariness in 



. . . , £ = £(<S). Then, ([5]) can be written as follows: 



d s $(S) ■ S + rf(S)9 + e- 



q(s) ■ ve + oz(s) = o. 



(8) 





d^(<Sii)+7?OSri) = 0; 

d v J(S u ) + e~ 1 q(Sn) = 0; (10) 
d a $(Sn) + t(S u ) = 0. 



Finally, when S = Sm, Sm fl Sm = {d,Vd}; this time, dH) is satisfied 
identically (if and only) if 



ip can only depend on a, 9, andVa; 



de^Sm) + ff(Sm) = 0; 

9dJ(Sm) + {dva$(Sm) + 9~ x q(S m )) ■ V9 + 9^(S m ) = 0. 



(11) 
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Note that the last relation in each of the three sets (|9])- (|lip may be thought 
of as defining the constitutive mapping that delivers the internal entropy 
production, in the following three forms: 

Cu(Su) = -dJn(S u ), (12) 
£hi(<Siii) = - dJ m (S m ) - 0' 1 (dv a $m(Sm) + ^ qm(Siii)) ■ V0, 

that is, in terms of the mappings delivering the heat influx and/or the free 
energyo 

In the appendices to [9] and |10j . Green and Naghdi assert that, when (|8|) 
is exploited to find restrictions on the constitutive mappings, the internal 
entropy production £ must be considered independent of the time rates of 
the state variables. This is indeed the case for their Type I theory, where 
Si n Si = 0, and they do arrive at the correct results Q in that case. But 
they regard that assertion about £ valid for Type II and Type III theories 
as well, whereas both Su n <Sn and 5m n «Sni are not empty; hence their line 
of reasoning is incorrect, as are their conclusions. 



Remark. It is instructive to compare the results (|9])- (|lljl with the results 
that are obtained if one makes about the internal entropy production, by 
far the most mysterious and least experimentally accessible of all the state 
functions, the minimally committing assumption that it is nonnegative: 

£ > 0, (13) 

an assumption that is at the conceptual core of standard continuum ther- 
modynamics. A direct consequence of (fl3j) is that (inequality ([5]) takes the 
familiar form ([6]) and) inequality (JHJ) reads: 

d s ^(S) ■ S + 7/(5)0 + 9" 1 q(S) • V0 < 0. (14) 

From (|14p . in the case of Type I theory one derives ([9|) 1 2 plus, in place of 
([IJlg, the classic inequality for the heat influx: 

q(Si) -V6<0. (15) 



4 Here, at variance with what we did when writing (p?))- (|lip to keep our notation lighter, 
we have carefully distinguished the heat-influx and free-energy constitutive mappings ac- 
cording to the type of theory they are meant for: e.g., we have written ipi in the case of 
Type I theory and '(/'ii, Vtii in the other two cases. 
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Likewise, for Type II theory one re-obtains (fTU|) 1 2 , while (fT0]) 3 is replaced 
by 

d a $(Sn) < 0; 

and, for Type III, the only change is that (jll j) 3 is replaced by the residual 
entropy inequality 

d a $(Sm) + e-^dvafcSm) + 6- 1 q(S m )) • V0 < 0. 

Green-Naghdi theories are given shape when the thermodynamic infor- 
mation in the lists from ([9]) to (|lip are inserted into the entropy balance 
dl]). In all cases, the same equation is arrived at if the same information 
is inserted into the energy balance ([3]), the procedure we adopt in the next 
three sections. 

3 Type I theory 

To begin with, we derive from definition ([!]), with the use of 

i) = and 77(0) = -d e ^(6), 

the following general expression for the internal energy: 

e = ij}(e)-edg$(e). 

Now, let us recall the classic prescriptions for free energy and heat influx: 

iP = -A0(log0 - 1) and q = -«V0, (16) 

where the heat capacity A and the thermal conductivity k are two positive 
constants (positivity of n is demanded for the Fourier law (|16p 2 to agree 
with the residual entropy inequality (fT5|) ). With (fT6]) 1 , one finds: 

r/ = A log and e = A0; (17) 

the end result is the well-known heat equation: 

X9 = KA9 + r, (18) 

At variance with (|16p . Green and Naghdi assume that 

M°) = ( 19 ) 
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and take 

£i(0,V0) = -£-V0, (20) 



for the entropy influx, 9q being a constant reference temperature; (|19p and 
([20]) imply that, respectively, 

W) = *T (21) 

and 

q l {Q,V0) = -K?-V6. (22) 
It follows from the above constitutive assumptions that 

= Ja^; (23) 
moreover, the first of (|12p specifies the internal entropy production: 

fi(0,v0) = ^-e^ve-ve. 

Insertion of (f23|) and ([22]) into (J3]) yields, with the use of Qg and ([2|) 2 , a 
balance equation pretty different from the classic heat equation (|18[) , namely, 



A0 = kA# + o (£i + s). 

Green and Naghdi get to (|18j) by assuming that ". . . the temperature 9 rep- 
resents departure from an equilibrium temperature 9q an d • • • time and space 
derivatives of 9 . . .are small of O(e) ..." ([TO], p. 197), which is slightly 
more than needed to take 

o fi(0,V0)~O, 6 s~r, 

and, said differently, amounts to restrict attention to almost isothermal and 
spatially uniform heat propagation processes, for which 

9 ~ 6 and |V0| ~ 0. (24) 

While one might be willing to study the nonlinear parabolic PDE Green 
and Naghdi arrive at before proceeding to their so-called 'linearization', 
namely, 

X9 = k (T x div (9 V0) + 9-% r, 
we see no point in deducing (I18p through the involved path they chose. 
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4 Type II theory 



We read in |10j the following recipe for Type II free energy: 

1 9^ 1 K* / 1 K* \ 

Ma, 9, Va) = --A — + ^ Va " Va ( = W ) + g ^ Va ' Va ) ' (25) 

where denotes a non-standard thermal conductivity coefficient. Note 
that, as a direct consequence of assuming that the free energy does not 
depend on the thermal displacement, it follows from (fTUjl o that the internal 
entropy production is null: 

&i(a,0,Va) = O@ (26) 

The entropy associated with prescription (|25|) for the free energy turns out 
to be the same as in ([2Tj) : 

rfn(a,e,Va) = X^- = rfi(e); (27) 

consequently, 

en(a,0,Va) = -A — + -— Va • Va. (28) 

A peculiar property of Type II theory is that the energy and entropy influxes 
are both determined by the free energy, which plays the role of a potential 
function; we have from (jlOP a that 

q n (a, 9, Va) = -—0 Va = 9 h u (a, 9, Va). (29) 

Remark. A well-known consequence of (jl6p 2 is that the Fourier influx vector 
q satisfies the classic heat- conduction inequality 

q- V0 < 0, 

which makes geometrically evident that heat influx and temperature gra- 
dient should always make an obtuse angle at any given body point in any 
given process of thermal propagation, in accord with common experience. 

5 In our opinion, this artifact should be considered less surprising than it seemed when 
[8] and [10] appeared, because of the number of ad hoc assumptions required to put it 
together. 
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One may ask whether the energy influx vector qn in (|29p would satisfy the 
same inequality. As a straightforward computation shows, 

q n {a, 9, Vq) • V9 = -^-9 ^|Va| 2 ; 

hence, provided that one assumes that k* > 0, in a Type II theory the 
energy influx and temperature gradient form an obtuse angle only when 
|Va| 2 grows bigger. 

Under the current assumptions, on making use of ()26|) . ()27|) . f|28|) . and 
([29|) 2 . both the entropy balance ([I]) and the energy balance ([3]) become: 

X9 = K*Aa + 9 s, 

On recalling that 9 = a, the previous equation looks as follows: 

Xa = K*Aa + 6»o s, (30) 

a hyperbolic PDE for a. Green and Naghdi do not follow this path in |10j : 
they rather differentiate (|30p with respect to time, presumably in order 
to get a PDE for temperature, a more familiar unknown field that thermal 
displacement, and a field for which assigning initial and boundary conditions 
is no problem. Following their suggestion, we get a hyperbolic PDE for 9 
that opens the way to study thermal wave propagation, namely, 

\0 = k*A6 + 6 s. (31) 
Unfortunately, the equation Green and Naghdi arrive at in |10j is: 

X9 = k* A9 + r,@ (32) 
an equation which coincides with (I3ip only if one takes 

r ~ 9q s, 

that is, only if attention is restricted to almost isothermal processes, as was 
done in the case of Type I theory by stipulating that (|24pi hold. 

6 Most probably because of a typographical error, the dot signaling time differentiation 
of the last term is missing in [10] . 
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5 Type III theory 



Among the three types of heat propagation theories Green and Naghdi pro- 
posed, this is in our opinion the one whose deduction is, so to speak, the 
most adventurous. Here is how we reconstruct it, to the extent that it is 
possible, from [SJ, the paper where Green-Naghdi present their Type III 
thermoelastic theory. 

The Type III free-energy mapping is the same as for a theory of Type 

II: 

e 2 k* 

V>m(a,0,Va) = -A — + — Va • Va. (33) 
As to energy and entropy influxes, it is proposed to take 

g m (a,0,Va, W) = - - - 9Va-^-9V9 = 9h m (a,9,Va,V9). (34) 

9o "o 

Note that 

q m (a, 9, Va, V0) = q Y (9, V0) + q„(a, 9, Va) - —9 Va, 

so that, in addition to k, and k*, another thermal conductivity modulus - 
but not another form of heat conduction! - enters the theory; presumably, 
k** > 0. In view of (|HP o. and as a consequence of (|33|) and (|34p . we have 
that 

Q 

%i(a, 9, Va) = A — , (35) 
9o 



so that, on recalling (|2Tj) and (|27j) . entropy, whatever the type of theory, 
depends only on temperature, and in a manner that can be regarded as a 
first approximation for 9 ~ 9q of the classic prescription (1171) x . Moreover, 
as a consequence of (j33j) and (f34"|) . we have that 

£n(a, 9, Va, V0) = 9~ l (^Va + £ w) • V0, (36) 
or rather, equivalently, 

6n(a, 6», Va, V6>) = 0" 1 (^-|(V«) 2 + £| W| 2 

Thus, once again, the expression for a quantity that should never be negative 
- this time, the internal entropy production - turns out to include a term 
that, just as the first one between brackets, could instead induce negativity 
during some process. 
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But this is not all. Green and Naghdi [9] claim that substitution of 
Eqs. (|35p . ()34P o and (j36[) into the 'linearized' form of what they refer to as 
'Eq. (2.17)' would yield: 

\e = KA6 + K*A9 + f. (37) 

Now, 'Eq. (2.17)' is not found in [S|. If, imitating the procedure followed 
in [SI [10] , we perform the same substitution in the entropy balance ([T]) , we 
find: 

X6 = kA8 + k*Aq + 9 s + k**Aq + o £ni(a, 0, Va, V6>), 
or rather, on accepting again the pseudolinearization (|24|) . 

A# ~ kA(9 + K*Aa + r + K**(Aa + Va • V0). 

We fail to see how one could get rid of the last two terms on the right side 
in a manner both convincing and different from taking k** = 0. But then, 
one could infer that, for Type III theory, the internal entropy production is 
approximately null: 

Sn(a, 9, Va, V0) ~ 0. (38) 

6 Conclusions 

In a number of papers appeared between 1977 and 1996 [U EJ EJ UHl 
[TT| [12], Green and Naghdi have proposed three types of heat propagation 
theory. According to one of these theories, heat propagates without internal 
entropy production; in two of them, heat waves are possible. For these and 
other reasons, all Green-Naghdi theories, and especially so Type III theory, 
have soon attracted an interest that does not seem to fade away even at 
the time of this writing, and have been applied in a number of disparate 
physical circumstances, where propagation of heat is coupled with elastic 
deformations of solids, flow of viscoelastic fluids, etc. (an updated account 
of the applications we allude at is found in [27]). Unfortunately, in our 
opinion, Green-Naghdi theories have the severe limitations pointed out in 
our present exposition; the three main ones are recapitulated hereafter. 

• The internal entropy dissipation £ is presumed to be the object of a 
constitutive prescription, which is in terms of the same state variables 
as Helmholtz free energy ip and energy and entropy influxes q, h, and 
turns out to depend on their constitutive prescriptions as specified by 
(fT2l) : in particular, £ must be null in all processes considered within the 
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framework of Type II theory. Now, a fundamental tenet of continuum 
thermodynamics, that was not made part of Green-Naghdi theories 
neither by their creators nor by their followers, is that £ must be non 
negative. To conform to this tenet for Type III theory, one should be 
able to prove that £ is non negative along whatever process, something 
that turns out to be possible under assumptions that would kind of 
emasculate that theory, at least partlylll 

• The energy and entropy influx vectors are taken proportional via the 
temperature (q = Oh), an assumption that has been proven to hold 
only when heat propagation is, in a precise sense, isotropic [T]. 

• The notion of thermal displacement plays a rather ambiguous role as 
a state variable, in that: (i) it does not enter the constitutive prescrip- 
tions per se, but only through its time and space derivatives a, Va 
and Vd; (ii) as a exemplified by (i30j) and (f32|) . the direct outcome of 
a theory, a PDE for a, is quickly differentiated with respect to time to 
arrive at a PDE for the temperature 9, a field variable for which is in- 
deed natural to lay down initial and boundary conditions (we wonder 
how one could be confident that an assignment of a at an initial time 
or on part of the boundary could be reproduced in the laboratory.) 

In spite of these criticisms, we believe that Green-Naghdi heat conduction 
theories, when amended as indicated, are worth studying and applying, es- 
pecially if a convincing statistical-mechanics interpretation is found for the 
thermal displacement. Thermal wave propagation is an example of worth- 
pursuing research theme: it would be interesting to compare the predictions, 
say, of equation (|37p with those of other mathematically similar equations 
that have been proposed (for an account of the literature up to 1990 see 
|144 I15j : see also [I], where other references are found, and a very recent 
contribution of one of us |25j). 
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